Gamow Teller strength in 54 Fe and 56 Fe 
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Through a sequence of large scale shell model calculations, 
total Gamow- Teller strengths (S+ and S-) in 54 Fe and 56 Fe 
are obtained. They reproduce the experimental values once 
the err operator is quenched by the standard factor of 0.77. 
Comparisons are made with recent Shell Model Monte Carlo 
calculations. Results are shown to depend critically on the 
interaction. From an analysis of the GT+ and GT— strength 
functions it is concluded that experimental evidence is con- 
sistent with the 3(N — Z) sum rule. 
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The charge exchange reactions (p, n) and (n,p) make it 
possible to observe, in principle, the total Gamow- Teller 
strength distribution in nuclei. The experimental infor- 
mation is particularly rich in 54 Fe and 56 Fe [Gl-p)) and 
the availability of both GT+ and GT— makes it possible 
to study in detail the problem of renormalization of err 
operators. Moreover, these nuclei are of particular as- 
trophysical interest ||, and they have been the object of 
numerous theoretical studies. 

In this paper we present the results obtained with the 
largest shell model diagonalizations presently possible. 
First we concentrate on the study of the total strengths 
S+ and S-. After a brief review of existing calculations, 
we estimate the exact values in a full Ohuj space, stressing 
the need of ensuring the correct monopole behaviour for 
the interaction. 

The second part of the paper deals with the GT+ and 
GT— stength functions. The analysis will confirm that 
the "standard" quenching factor of 0.77 is associated to 
suppression of strength in the Ofuv model space, but that 
little strength is actually "missing" (i.e., unobserved). 

I. Total Strengths S + and 5 

The experimental situation is the following: 

i) 54 Fe(n,p) 54 Mn, 

S+ = 3.1±0.6 from @, strength below 10 MeV. 
S + = 3.5±0.3±0.4 from g], strength below 9 MeV. 

ii) 54 Fe(p, n) 54 Co, 

S- = 7.5±1.2 from fjj, strength below 15 MeV. 
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S- = 7.8±1.9 from [3], strength below 13.5 MeV. 
S- = 7.5±0.7 from gj, strength below « 24 MeV, 
(but see discussion in section II). 

iii) 56 Fe(?i,p) 56 Mn, 

S+ = 2.3±0.2±0.4 from j|, strength below 7 MeV. 
S+ = 2.9±0.3 from g, strength below 8.5 MeV. 

iv) 56 Fe(p, n) 56 Co, 

S- = 9.9±2.4 from @, strength below 15 MeV. 

The theoretical approaches include: shell model calcu- 
lations in the pf shell with different levels of truncation, 
RPA, quasiparticle RPA and Shell Model Monte Carlo 
(SMMC) extrapolations. Let us examine the results. 

54 Fe. Previous shell model calculations. Throughout 
the paper / stands for I/7/2 and r for any of the remain- 
ing orbits. Truncation level t means that a maximum of t 
particles are promoted from the I/7/2 orbit to the higher 
ones, i.e., that the calculation includes the following con- 
figurations: 

tn— noj.no jn—n Q -l r n + l j?n—na—ty,no+t 

with n = A — 40. no is different from zero when more 
than eight neutrons (or protons) are present and at t = 
n — no we have the full space calculation. 

Given a choice of t — t p for a parent state having no = 
nop, to ensure respect of the 5_ — 5+ = 3 (AT — Z) sum 
rule, the truncation level for daughter states having no = 
nod, must be taken to be td — t p + 1 + no p — nod- 

In the simplest case we have t — i.e., 0p-2h con- 
figurations with respect to the 56 Ni closed shell for 
the 54 Fe ground state, and lp-3h configurations for the 
54 Mn daughters. The result —S+ =10.29 Gamow Teller 
units — is independent of the interaction. The calcula- 
tion was extended to t — 1 by Bloom and Fuller j7| , using 
the interaction of ref. H, obtaining S+ — 9.12. A similar 
calculation by Aufderheide et a l. S | yields S+ — 9.31 (in- 
teraction from JlO| , p| ). Muto p2fj using the interaction 
p3| , made a t =2-like calculation that did not respect 
the 3 (AT — Z) sum rule, but the author estimated the in- 
fluence of this violation and proposed S+ =7.4. Finally 
Auerbach et al. E3 have made a t =2 calculation using 
the interaction, MSOBEP, fitted in |J| (BR from now 
on), and obtain S+ =7.05. 

QRPA. The calculation of Engel et al. |§ yields S+ = 
5.03, to be compared with QRPA or RPA calculations of 
Auerbach @ leading to 5+ = 6.70. 

Shell Model Monte Carlo. The calculation of Alhassid 
et al. |]l7j , in the full pf shell, using the BR interaction 
extrapolates to S+ = 4.32 ± 0.24. (here the error bar 
includes only the statistical uncertainties but not those 
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associated to the extrapolation or to possible sistematic 
errors of the method) . 

Large t Shell Model calculations. All the previous re- 
sults point to a reduction of GT+ strength as correla- 
tions are introduced and to a rather large dispersion of 
the calculated values depending on the interaction and 
the approach used. Therefore, to obtain a reliable S+ 
value, the method and the interaction must demonstrate 
their ability to cope with a large number of other prop- 
erties of the region under study. Calculations in the pf 
shell 1 18 1 using the KB3 interaction — a minimally mod- 
ified version of the Kuo Brown G-matrix — fulfill 
this condition since they give an excellent description of 
most of the observables in the region up to A=50. The 
same interaction was used years ago in perturbation the- 
ory to describe nuclei up to 56 Ni ^0|, with fair success. 
It should be mentioned that the monopole modifications 
in KB3 involve only the centroids Vfj and V/ r . The V rr 
values were left untouched and may need similar changes. 

It is not yet possible to perform a full pf shell calcu- 
lation in 54 Fe. However, we can come fairly close by fol- 
lowing the evolution of the total strength as the valence 
space is increased. The shell model matrices are built and 
diagonalized and the GT strengths calculated with the 
code ANTOINE Full Lanczos iterations in spaces 

that reach maximum m-scheme dimension of 1.4 ■ 10 7 are 
necessary for the parent states. Acting on them with the 
err operator to calculate the strength, leads to spaces of 
m-scheme dimension of 4.1 • 10 7 . 

In addition to KB3, to compare with the results of 
the SMMC extrapolations of |17|, we have used the BR 
interaction |Tq| . The results are collected in table I and 
we proceed to comment on them. 

1. The t =5 calculation should approximate the ex- 
act ground state energy reasonably well, as can be 
gathered from the small gain of 270 KeV achieved 
when increasing the space from t =4 to t =5. 

2. The SMMC result using the BR interaction, 
—55.5 ± 0.5 MeV, is some 1 MeV above the exact 
energy since our t =5 result gives an upper bound. 
Consequently the SMMC error bars in |17| are un- 
derestimated. 

3. The result of our t =2 calculation using BR differs 
slightly from the one in Q (7.22 vs. 7.05). This 
is due to the readjustment of the single particle 
energies made in |1J] with respect to the values of 
©■ 

4. Auerbach et al. proposed an extrapolation of 
their t = 2 calculation to the full space, based 
on the behaviour of S+ in 26 Mg as a function of 
B(E2; 2+ -> 0+). Although it is true that there is 
a qualitative correlation between these two observ- 
ables (the bigger the quadrupole collectivity the 
smaller the S+ value), it is difficult to go further 
and to obtain a quantitative prediction. In figure 1 



we have plotted the S+ values vs. B(E2) for the BR 
interaction and several truncations. It is clear that 
no simple correlation pattern comes out. Notice 
also that the extrapolation in gives 5+ = 6.4 
compared to 5*+ = 5.5 in the t — 5 calculation. 

Before we discuss the differences between the results 
of the BR and KB3 interactions and between shell model 
diagonalizations and Monte Carlo extrapolations we ex- 
amine the situation in 56 Fe. 

56 Fe. Bloom and Fuller made a t — calculation 
that yields S+ = 10.0 (interaction from ||). Anantara- 
man et al. [Q (interaction from |U^,[nj) obtain 5+ = 
9.25 for t = and S+ = 7.38 for t = 1. The SMMC 
result ]23[ | is shown in table I together with the numbers 
coming out of several truncations for both KB3 and BR 
interactions. 

The influence of the interaction. The interactions 
KB3 and BR lead to different single particle spectra for 
57 Ni. The sequence of levels obtained in the calcula- 
tions up to t — 4 are compared with the experimen- 
tal data in table II. It is apparent from the table that 
the BR interaction places the I/5/2 orbit too low. As 
a consequence the dominant configuration in the 56 Fe 
ground state predicted by BR is (l/V/2) 14 (l/5/2) 2 instead 
of (I/7/2) 14 (2^3/2) 2 as given by KB3. This explains the 
very large difference in S+ values observed in table I, 
already at the t =0 level: For a pure (l/V/2) 14 (2p3/2) 2 
configuration the total strength is 10.3, while for a pure 
(l/7/2) 14 (l/5/2) 2 h amounts to only 5.7. In 54 Fc the 
situation is not so dramatic because the leading configu- 
ration, (I/7/2) 14 , is the same in both cases. Still the BR 
value is 20% smaller than the KB3 one, due to an excess 
°f I/7/2 - I/5/2 mixing in the ground state. From that 
we conclude that the BR interaction underestimates the 
S+ values for nuclei with N or Z greater than 28. 

Table II also shows the valus of the "gaps" defined by: 

A = 2B£( 56 Ni) - BE( 57 m,l~) - BE{ 55 Wi). 

The strong staggering between even and odd values 
of t makes it difficult to obtain a reliable extrapolation. 
The overall trend for the gap is to decrease as t increases. 
Nevertheless, it is probable that the exact value for KB3 
will remain somewhat larger than the experimental one. 
In this case a slight revision of the monopole terms would 
be needed. 

Shell Model extrapolations. In figure 1 we show 
the evolution of the total GT strength with the level of 
truncation in 54 Fe , 56 Fe and several cases ( 48 Ti, 50 Ti, 
48 Cr, and 50 Cr) for which exact results are available. If 
we continue the 54 Fe calculated values with lines paralel 
to the A=50 ones, we get the following extrapolated val- 
ues: 

54 Fe ; S+(KB3)= 6.0 ; S+(BR)= 5.0 
If we assume that the value of the difference between the 
t =4 to t =5 result and the exact one is the same in 54 Fe 
and 56 Fe, the corresponding extrapolation is 

56 Fe ; S+(KB3)= 4.5 
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these values are fully consistent with the experimental 
results if we use the standard 0.77 renormalization of the 
Gamow- Teller operator ( p^-p6| and section II) . For 54 Fe 
we have 

S+(exp)= 3.1±0.6 ; 3.5±0.7 vs. 

S+(KB3)= 3.56 ; 5 + (BR)= 2.96 
the corresponding predictions for S- are compatible - 
again within the 0.77 renormalization — with the exper- 
imental results 

SL(exp)= 7.5±1.2 ; 7.8±1.9 ;7.5±0.7 vs. 

SL(KB3)= 7.11 ; S_(BR) = 6.52. 
If we turn to 56 Fe the corresponding numbers are: 

S+(exp)= 2.3±0.6 ; 2.9±0.3 vs. S"+(KB3)= 2.7, and 

S_(exp)= 9.9±2.4 vs. S_(KB3)= 9.8. 
We have prefered to omit error bars in the extrapolated 
values which are simply reasonable visual guesses. How- 
ever, they fall so confortably in the middle of the ex- 
perimental intervals, that an estimate of computational 
uncertainties would leave the conclusions unchanged. 

Comparison of Monte Carlo and Shell Model 
extrapolations. The differences between our results 
and those of ref. |^3| are mostly — but not only — due 
to the use of different forces. With the same force the 
shell model extrapolations yield values that are some 20% 
larger than the SMMC ones. The discrepancy is probably 
related to the lack of convergence of the SMMC energies 
detected in table I. 

Note J27|. In the most recent SMMC calculations 
with finer Af3 steps of 1/32 (instead of 1/16) the bind- 
ing energy goes down by 1 MeV thus eliminating the 
problem mentioned in point 2 above. Furthermore S+ 
becomes 4.7±0.3 in full agreement with our extrapolated 
value . SMMC values have also become available for the 
KB3 interaction [§f]. The values for 54 Fe and 56 Fe are 
6.05±0.45 and 3.99±0.27, again in very good agreement 
with our values. 

II. Strength functions: standard quenching and 
missing intensity 

In fig. 3 we show the total I = cross sections obtained 
by Anderson et al. for 54 Fe(p, n) 54 Co. The individual 
peaks in table I of |J] have been associated to gaussians 
of cr=87 keV (the instrumental width) for the lowest and 
cr=141 keV for the others. The "background" (the area 
under the dashed lines) is obtained by converting the 2 
MeV bins in table III of j| into gaussians of o=1.41 MeV. 
Since there is no direct experimental evidence to decide 
how much of this background is genuine strength, two 
extreme choices are possible to extract S- : either keep 
the whole area in the figure (i.e., S- = 10.3±1.4), or only 
the area over the dashed line (i.e., S- = 6.0 ± 0.4). An 
intermediate alternative consists in keeping what is left 
of the background after subtracting from it a calculated 
contribution to quasi free scattering (QFS). The result- 
ing profile (with S- = 7.5 ± 0.7, the number adopted in 
I.ii) is shown in fig. 4 (tables I+II of B) and compared 
with the Lanczos strength function (see [£6| for instance) 



obtained after 60 iterations in a t = 3 calculation for the 
parent state and t — 4 for the daughters (the peaks are 
broadened by gaussians of cr=87 keV for the lowest, and 
cr=212 keV for the others). The areas under the mea- 
sured and calculated curves are taken to be the same (we 
know that upon extrapolation to the exact results they 
coincide). To within an overall shift of some 2 MeV, the 
two profiles agree quite nicely. The discrepancy is easily 
traced to the (too large) value of the gap in table II at 
this level of truncation. 

Although a calculation closer to the exact one would 
be welcome, the elements we have point to a situation in 
all respects similar to that of the 48 Ca(p, n) 48 Sc reaction, 
analyzed in |26|. What was shown in this reference can 
be summed up as follows: 

• The effective trr operator to be used in a Qhu calcu- 
lation is quenched by a factor close to the standard 
one (« 0.77) through a model independent mecha- 
nism associated to nuclear correlations. 

• The rest of the strength must be carried by "intrud- 
ers" (i.e., non 0hu> excitations). Only a fraction of 
this strength is located under the resonance, but 
intruders are conspicuously present in this region 
and make their presence felt through mixing that 
"dilutes" the OHuj peaks causing apparent "back- 
ground" . 

In all probability, the long tail in fig. 3 corresponds 
to intruder strength and should be counted as such. 
What is achieved by subtracting the QFS contribution 
amounts — accidentally but conveniently — to isolate the 
0hu> quenched strength. It is to this contribution that 
the notion of standard quenching applies but it should 
be kept in mind that the remaining strength — necessary 
to satisfy the 3 (TV — Z) sum rule — is not missing, but 
most probably present in the satellite structure beyond 
the resonance region as hinted in the very careful analysis 
of ref. §j. ' 

In fig. 5, to make a meaningful comparison with the 
54 Fe(n,p) 54 Mn data of Q the spikes of a t = 3 calcu- 
lation have been replaced by gaussians with a = 1.77 
MeV, chosen to locate some strength at —2 MeV, where 
the first experimental point is found. The resulting dis- 
tribution is then transformed into a histogram with 1 
MeV bins. The agreement is quite satisfactory. It should 
be pointed out that the measures of [Q are displaced to 
lower energies by some 700 keV with respect to those of 
P). Otherwise, the experiments are in good agreement, 
and both show satellite structure beyond the resonance 
(not included in fig. 5, but visible in figs. 10 and 7 in jl]] 
and respectively. 

Finally, in figs. 6 and 7 we show the corresponding 
results for 56 Fe targets, for which a t = 2 truncation 
level was chosen, going to t = 4 for (p, n) , and t = 2 for 
{n,p). Though this numerical limitation is rather severe, 
the agreement with the data remains good enough to 
support the main conclusion of this paper: 
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The GT strength functions for 54 Fe and 56 Fe, in the 
resonance region and below, are well described by QTilj 
calculations that account for (0.77) 2 of the total strength. 
The remainder, due to intruder states, is likely to be 
present in the observed satellite structures, so that the 
3(N — Z) sum rule is satisfied. 

We have also shown that spurious reductions of the 
GT+ strength can occur due to defects of the effective 
interaction as it is most probably the case for some of the 
results of ref . pj| . 

This work has been partly supported by the IN2P3 
(France) - CICYT (Spain) agreements and by DGI- 
CYT(Spain) grant PB93-264. 
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TABLE I. 54 Fe ^ 54 Mn and 56 Fe ^ 56 Mn Gamow Teller 
strength S+ in units of the GT sum rule. In column 4 nuclear 
two body energy of the ground state of 54 Fe (in MeV). 



KB3 3/2 5/2 1/2 A 


BR 3/2 5/2 1/2 A 


t = 0.0 0.38 1.15 8.57 
t = 1 0.0 0.47 1.14 7.33 
t = 2 0.0 0.72 1.16 8.10 
t = 3 0.0 0.76 1.14 7.74 
t = 4 0.0 0.86 1.14 7.90 


0.48 0.00 3.06 7.42 
0.07 0.00 2.11 5.80 
0.07 0.00 2.27 7.01 
0.00 0.08 1.89 6.41 
0.00 0.11 1.83 7.21 


EXP 0.0 0.77 1.11 6.39 





TABLE II. Excitation energies of the low-lying states in 
57 Ni and the gap A in MeV (see text). KB3 and BR results for 
several truncations, compared with the experimental results. 
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FIG. 1. 54 Fe. S+ (in units of the GT sum rule) vs. B(E2) 
(2+ -> 0+) (in e 2 /m 4 ). BR interaction. The last point to the 
right is extracted from the results of |17{. 
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FIG. 3. Experimental GT- strength in 54 Fe 
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FIG. 2. Total GT+ strength as a function of the truncation 
level t. Calculated values are linked by lines, the isolated ones 
are extrapolations. KB3 interaction, in all cases but one. 
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FIG. 4. Experimental GT- strength in 54 Fe after QFS sub- 
traction |4j, compared toi = 3— >t = 4 calculations. 
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FIG. 5. Experimental GT+ strength in 54 Fe compared 
to t = 3 — * t = 3 calculations. 
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FIG. 7. Experimental GT+ strength in 56 Fe HU .compared 
to t = 2 — * t — 2 calculations. 
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FIG. 6. Experimental GT— strength in 56 Fe j^, compared 
to t = 2 — > i = 4 calculations. The IAS peak has not been 
removed from the data, but it is not included in the calcula- 
tions. 
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